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The <p meson productions in Au + Au and/or Pb + Pb collisions at AGS, SPS, RHIC, and LHC 
energies have been studied systematically with a hadron and string cascade model LUCIAE. After 

(N 

considering the energy dependence of the model parameter a in string fragmentation function 



and adjusting it to the experimental data of charged multiplicity to a certain extent, the model 
predictions for (j> meson yield, rapidity, and/or transverse mass distributions are compatible with 
the experimental data at AGS, SPS and RHIC energies. A calculation for Pb + Pb collisions at 
LHC energy is given as well. The obtained fractional variable in string fragmentation function 
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shows a saturation in energy dependence. It is discussed that the saturation of fractional variable 
in string fragmentation function might be a qualitative representation of the energy dependence of 
nuclear transparency. 
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Strangeness enhancement was suggested in the early eighties [|l| as one of the most promis- 
ing signatures for the creation of a Quark-Gluon Plasma (QGP) phase in relativistic nuclear 
collisions. Following the experimental observations on strangeness enhancement in proton- 
nucleus and nucleus-nucleus collisions at the SPS energies 0, |3], f|, ^ ^ [7J the WA97 has 
measured a clear enhancement of multi-strange baryons (A, H, O) with their strange quark 
content in 158A GeV/c Pb + Pb collisions relative to p + Pb collisions ||. Recently the 
STAR data on the strangeness production in Au + Au collisions at v /s nn =130 GeV were 
reported ||. 

TABLE I: Global hadron multiplicity in relativistic nucleus-nucleus collisions 



Reaction 


Au + Au 


Pb + Pb 


Au + Au 


Pb + Pb 


y/s^, (GeV) 


4.66 


17.3 


130 


5500 


Centrality 


Exp. 


<4% 


< 4% 


< 5% (BRAHMS) 


< 10% 


LUCIAE 


b <2.82 fm 


b <3.5 fm 


b <3.2 fm 


b <4.54 fm 


N ch 


Exp. 
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4114 
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(< 7T+ > + < 7T _ >)/2 


Exp. 


105 3 ) 


611 4 ) 






LUCIAE 


102 


611 




5851 2 ) 


a 


0.4 


1.3 


12 


18 


< z > 


0.286 


0.565 


0.923 


0.947 



1. -4.7< r] <4.7, from (2£ 
3. 0.6 < y < 2.6, taken from pC 



2. in full phase space 
4. cf. 



As the mesonic counterpart, the enhancement of (j) meson production in relativistic nuclear 
collisions was also suggested as an evidence of the QGP formation in Ref. |i~0|| , since in the envi- 
ronment of a QGP the copious strange and antistrange quarks originating from gluon annihilation 
would be very likely to coalesce forming <p mesons during the hadronization period. Due to the 
small cross sections of 4> mesons interacting with non-strange hadrons EJ, [l0|], penetrating cj) mesons 
are also messengers of the early stage of the colliding system. Thus, the <fi meson is not only a 
promising signature for the QGP formation but also a good probe to study the reaction dynamics. 



After the experimental observations on <ft meson productions in Au + Au collisions at AGS [11] 



and sulfur- nucleus collisions and Pb + Pb collisions at SPS [12, 13, 14 1, the STAR collaboration 

2 



11 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

(158A GeV/c)Pb+Pb 



5 
Z 



NA49 

NA49 (mirror) 
LUCIAE 

i i i i I i i i i I i i 



> 

a 

T3 

■a 

s" 



z 

'a 



10 



10' 



10 



10"' 







fA (158AGeV/c)Pb+Pb : 


- a, p.<y<j.S) 


- 


A • 


- 


I A 

• 




A 




A • 




" A i 




A 




• NA49 




a LUCIAE 




mmImmI 





0.5 1 

m t" m o ' GeV ^ 



FIG. 1: Rapidity (left panel) and transverse mass (right panel) distributions of 4> mesons in Pb+Pb 
collisions at 158 A GeV/c 
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FIG. 2: Transverse mass distributions (-0.5< y <0.5) of (f> mesons in Au+Au collisions at ^/s nn =130 
GeV 



The model studies on the d> meson enhancement in relativistic nucleus-nucleus collisions are rare 



1 16, 17, IS, 19 1 and there exists up to now no theoretical investigation on the energy dependence of 
<fi production from AGS, to SPS, to RHIC, and up to LHC energies systematically. In this letter, a 
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hadron and string cascade model, LUCIAE [£0| was employed to investigate first time the energy 
dependence of <\> meson production systematically. We have successfully used LUCIAE to study 
the enhanced production of multi-strange baryons (A, H, 0) and determined the model parameters 
related to the production of strange particles [21, 22]. 

TABLE II: <f> mesons yield in relativistic nucleus-nucleus collisions 



Reaction 


V ( s nn) 

(GeV) 


Centrality 


<N^> 


Exp. 


LUCIAE 


Exp. data 


LUCIAE 


Au + Au 


4.66 


301 1 ) 


307 1 ) 


0.252±0.107 2 ) 


0.14 


Pb + Pb 


17.3 


<4% 


b <3.5 fm 3 ) 


7.6±1.1 3 ) 


6.48 


Au + Au 


130 


< 11% (STAR) 


b <4.31 fm 


5.73±1.06 4 ) 


5.10 


Pb + Pb 


5500 


< 10% 


b <4.54 fm 




132 5 ) 



1. < Np art > 

3. cf. 0] 

5. full phase space 



2. within 0.9< y <1.4 |TT] 
4. ^ within -0.5< y <0.5 



The LUCIAE model is based on FRITIOF 23], which is an incoherent hadron multiple scatter- 
ing and string fragmentation model. In FRITIOF, the nucleus- nucleus collision is depicted simply 
as a superposition of nucleon-nucleon collisions. What characterizes LUCIAE beyond FRITIOF 
are the following features: First of all, the rescattering among the participant and spectator nu- 



cleons and the produced particles from string fragmentation are generally taken into account [24]. 
However, as proposed in [|l], [T(| the effects of the final state interactions on the (j) meson production 
and propagation are neglected. Secondly, the collective effect in the gluon emission of strings is 
considered by firecracker model |2q| . Thirdly, a phenomenological mechanism for the reduction 
of the s quark suppression in the string fragmentation is introduced |2l| resulting the effective 
string tension and therefore the pertained JETSET parameters vary automatically with the en- 
ergy, centrality, and size of reaction system. Those JETSET parameters: parj(2), (3), (1), and 
(21) in program are, respectively, the suppression of s quark pair production in string fragmenta- 
tion compared with u or d pair production, the extra suppression of strange diquark production 
compared with the normal suppression of strange quarks, the suppression of diquark-antidiquark 
pair production compared with quark-antiquark pair production, and the width in the Gaussian 
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FIG. 3: (a) fractional variable in string fragmentation function < z > and (b) < Nqh >/\/ s AA as 
a function of yjs nn 

From the AGS energy ~ 5 GeV) to SPS ( v /s^ ~ 20 GeV) to RHIC ( v /s^ ~ 200 

GeV) and then to the LHC energy (-y/ s nn ~ 5500 GeV) there are four energy magnitudes spanned. 
Among the nucleus-nucleus collisions in above energy region their dynamic behaviors, of course, 
must be different observably from each other. To describe the nucleus-nucleus collisions in such 
a wide energy region one either selects suitable model for different energy or considers the energy 
dependence of model ingredients properly if one single model is used. This idea has indeed been 
adopted in e + e~ studies quite early [[^j. Recently the thermal model with two main parameters 
(temperature T and baryon chemical potential was used to fit the experimental data of global 
hadronic yields for nucleus-nucleus collisions at SIS (y/s nn ~ 1 GeV), AGS, SPS, and RHIC en- 



ergies, the obtained values of T and fis are quite different from each other [28]. Therefore, for 
dynamical models it is also needed to consider the energy dependence of the involved dynami- 
cal ingredients properly in order to describe well the nucleus-nucleus collisions in region spanned 
four energy magnitudes. Of course, some dynamical ingredients might be sensitive to the energy 
and others might not be. Many dynamical models (including LUCIAE) working well for nucleus- 
nucleus collisions at SPS energy, however, overestimate the charged multiplicity at RHIC energy. 
One reason might just be that there are ingredients either energy independent or energy dependent 
without taken into account properly. 

In LUCIAE model (in FRITIOF and in JETSET eventually) a very important dynamical in- 
gredient, the string fragmentation function, is energy independent. Thus the LUCIAE model is 
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extended in this letter to take the energy dependence of string fragmentation function into account. 



As mentioned in [26] that the string fragmentation function f(z), which expresses the probability 
that a given fractional variable z is sampled, could be arbitrary in principle. In f26|| several such 
choices are given, besides the one of Lund fragmentation function. For describing easily the rel- 
ativistic nucleus-nucleus collisions over four energy magnitudes one of the simplest option for the 



string fragmentation function in 26] is employed here for the moment. That string fragmentation 
function reads 

/(*) = "(v 7 ^) x z<^~\ (1) 

where the fractional variable z refers to the fraction of light-cone momentum taken away by pro- 
duced hadron from the fragmenting string, the parameter a in this string fragmentation function 
is assumed to be energy dependent now. 

We first ran LUCIAE for Pb+Pb and Au+Au collisions, respectively, at SPS and RHIC energies. 
In LUCIAE calculations the centrality cut and the windows of rapidity (pseudorapidity) and pt 
were the same as those in corresponding experiments (the same later). Adjusting the parameter 
a(y/s nn ) in order that the global hadron multiplicity from LUCIAE could be compatible with the 
corresponding experimental data, the obtained LUCIAE results are given in Tab. 1 together with 



the experimental data. The data of Pb + Pb collisions at SPS energy were taken from [14] and 
Au + Au collisions at ^/s nn =130 GeV from p9|] . Since in LUCIAE the reduction mechanism of s 
quark suppression was introduced under the requirement that the pertained JETSET parameters 
mentioned above, must be reduced to the corresponding default values for the p+p collisions at SPS 
energy, the LUCIAE model with this mechanism works for SPS energy and above only. Therefore 
we ran LUCIAE without reduction mechanism of s quark suppression for Au + Au collisions at 
AGS energy but a empirical value of s quark suppression factor (parj(2)=0.45) extracted from 
experiments 28] was used in stead of default ones in JETSET. The experimental data for Au + Au 



collisions at AGS energy were taken from |3C]. The LUCIAE calculation for Pb + Pb collisions at 
LHC energy was also given in Tab. 1, in that calculation the a(y/s nn )=18, larger properly than 
that in Au + Au collisions at RHIC, was assumed. In Tab. 1 the fixed a(y/s nn ) parameter and the 
obtained fractional variable in string fragmentation function < z > (averaged over events), were 
given as well. 

Then we employed the parameter a{^s nn ) fixed above to calculate the <p meson yield, rapidity, 
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and/or transverse mass distributions in relativistic nucleus- nucleus collisions. The LUCIAE results 
of <p meson yield were compared with the corresponding data in Tab. 2, where the experimental 



data were taken from flll|| , ||T4f| , and [15] for Au + Au collisions at AGS energy, Pb + Pb at SPS, and 
Au + Au at ^/s n „=130 GeV, respectively. One sees in Tab. 2 that globally speaking the LUCIAE 
results are well compatible with experimental data. 

In Fig. 1 we compare the LUCIAE results of 4> meson rapidity (left panel) and transverse mass 
(right panel, 3.0<y<3.8) distributions with the experimental data in Pb + Pb collisions at SPS 
energy (data taken from |14|). It is shown in Fig. 1 that the LUCIAE model works somewhat 
better for the rapidity distribution than the transverse mass distribution in Pb + Pb collisions at 
SPS energy. The comparison between LUCIAE results and experimental data of <j) meson transverse 
mass distributions (-0.5<y<0.5) in Au + Au collisions at ^/s nn =130 GeV was given in Fig. 2 where 



data were taken from [15]. One sees in Fig. 1 and 2 that for transverse mass distribution LUCIAE 



model works better at RHIC energy than SPS since more hard 4> mesons are produced at RHIC 



than SPS due to the Schwinger mechanism |H]] in particle production of string fragmentation [18]. 

Fig. 3 (a) gives the fractional variable in string fragmentation function, < z >, as a function 
of \/s nn from LUCIAE calculations above. A saturation structure at around RHIC energy is 
observed in Fig. 3 (a) . Since the fractional variable z is the fraction of light-cone momentum taken 
away from the fragmenting string by the produced hadron, it must relate strongly to the nuclear 
transparency in nucleus- nucleus collisions. The formation time in which string is excited after 
hard scattering and the hadronic rescattering were considered as dominative factors in baryon 
stopping (nuclear transparency ) |32]] . However, one knows that the higher reaction energy the 
longer formation time, and then the less rescattering and higher nuclear transparency. As for 
the rescattering itself, the higher reaction energy leads to the more produced particles, and then 
the more rescattering and lower nuclear transparency. Thus the effects of formation time and 
the hadronic rescattering on the nuclear transparency cancel each other in certain extent. One 
might expect that the fractional variable in string fragmentation function is dominant in nuclear 
transparency and its saturated energy dependence is a qualitative representation of the energy 
dependence of nuclear transparency. 

The corresponding LUCIAE results of charged multiplicity per unit reaction energy 
< N c h > j ' yJsAA in full phase space as a function of ^s nn are given in Fig. 3 (b). A trend of 
saturation seems to be there in Fig. 3 (b) as well. One sees in Fig. 3 (a) that the hadron is 



created with higher fractional variable as y/s nn is increased and the fractional variable approaches 
almost to its maximum value one when yjs nn gets the LHC energy. On the other hand, Fig. 3 (b) 
denotes that the efficiency rate of energy to hadron (matter) is decreasing with increase of yjs nn 




and tends also to some kind of saturation. Since the fractional variable equal to one means that 
only one hadron would be created from fragmentation of an exited string and the number of excited 
strings (nucleon-nucleon collisions) is limited, the fractional variable approaching to unit must lead 
to the saturated trend in < 7V ch > / s/saa- Those saturation trends might be an indication of the 
saturation in particle production. 

In summary, we employed a hadron and string cascade model, LUCIAE, to study the (ft meson 
production in Au + Au and/or Pb + Pb collisions at AGS, SPS, RHIC, and LHC energies. By 
considering the energy dependence of the model parameter a in the string fragmentation function 
and adjusting it to the experimental data of charged multiplicity to a certain extent, the model 
results of (ft meson yield, rapidity, and/or transverse mass distributions were compatible with the 
corresponding experimental data for AGS, SPS and RHIC energies. A calculation for the (ft meson 
production in Pb + Pb collisions at LHC energy is also given. The obtained fractional variable in 
string fragmentation function shows a saturation in energy dependence at around RHIC energy. It 
is discussed that this energy saturation phenomenon could be a qualitative representation of the 
energy dependence of nuclear transparency. 
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